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Collider searches for dark matter in events with a Z boson and missing energy 
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Searches for dark matter at colliders typically involve signatures with energetic initial-state radi- 
ation without visible recoil particles. Searches for mono-jet or mono-photon signatures have yielded 
powerful constraints on dark matter interactions with Standard Model particles. We extend this 
to the mono-Z signature and reinterpret an ATLAS analysis of events with a Z boson and missing 
transverse momentum to derive constraints on dark matter interaction mass scale and nucleon cross 
sections in the context of effective field theories describing dark matter which interacts via heavy 
mediator particles with quarks or weak bosons. 
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The particle nature of dark matter is one of the 
greatest outstanding mysteries of cosmology and parti- 
cle physics. A suite of dedicated experiments seek to 
shed light on how dark matter interacts with Standard 
Model (SM) particles by looking for detection of ambient 
dark matter either directly scattering with heavy nuclei 
or indirectly through its annihilation into high energy 
Standard Model particles. An important third pillar to 
the search for the particle nature of dark matter is fur- 
nished by high energy particle accelerators, which can 
produce pairs of dark matter particles, which are ex- 
pected to manifest as an excess of events showing an im- 
balance in momentum conservation. Searches for missing 
transverse momentum are a major activity at the LHC 
precisely because of their potential connection to dark 
matter [l]. 

Searches for dark matter in missing momentum chan- 
nels can be classified based on the visible particles against 
which the invisible particles recoil. Existing experimental 
studies have considered cases in which the visible radia- 
tion is a jet of hadrons (initiated by a quark or gluon) |2l- 
H], a photon [HUlj, or a boson decaying into leptons [7]. 
These studies are performed in the context of an effec- 
tive field theory (EFT) which captures the physics of 
a heavy particle mediating an interaction between dark 
matter and quarks and/or gluons. In this article, we ex- 
tend the menu of such searches to include the case of a Z 
boson decaying into a pair of charged leptons (electrons 
or muons) and recast the recent ATLAS measurement 
of ZZ — 7> Mi/v [S] into a bound on production of dark 
matter in association with a Z boson^. Since our sig- 
nature consists of a pair of leptons consistent with a Z 
boson decay recoiling against transverse momentum car- 



This signature has previously been considered in the context of a 
coUider search for Z' decaying to invisible modes (Ref. ^Sj), and 
more recently in the dark matter context with a slightly more 
model-dependent framework in Ref. 1101 . 



ried by particles invisible to the detector, we refer to our 
selection as a "mono-Z" signature. 

We work in the context of EFTs where the dark mat- 
ter's primary interactions are with quarks or directly with 
electroweak bosons. In the case where interactions are 
primarily with quarks, the mono-Z signature arises from 
a Z boson which is radiated from a qq initial state, much 
like mono-jets or mono-photons. Such interactions also 
imply (depending on the specific form of the interaction) 
large rates for scattering with heavy nuclei. The case of 
direct interactions with a pair of Z bosons is more chal- 
lenging to connect to direct detection (however, see [TT|). 
and is a particular strength of collider searches. As usual, 
in both cases the EFT will break down at an energy not 
far above the one which characterizes the strength of the 
interactions written as higher dimensional operators, and 
is only a good description of the physics for processes tak- 
ing place at energies well below this cut-off scale. 

EFFECTIVE FIELD THEORY 
Interactions with Quarks 

Effective field theories for dark matter interacting pri- 
marily with SM quarks have been considered in Refs. |12| - 
m]. We consider the interactions, 

^ ( mn _ _ 1-,, 
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where x is the dark matter particle, which we assume to 
be a Dirac fermion, g is a SM quark, and the coefficients 
A parameterize the coupling strength of scalar (Dl), vec- 
tor (D5), axial- vector (D8), and tensor (D9) interactions 
between the two. The labeling scheme is adopted from 
Ref. [TH], with the choices dictated as those operators 
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which lead to a non-vanishing scattering rate with nu- 
cleons at small momentum transfer. We will typically 
consider one interaction type to dominate at a time, and 
will thus keep one A finite while the rest are sent to in- 
finity and decoupled. These operators are normalized so 
as to be consistent with minimal flavor violation. 




Interactions with Z Bosons 



One may also construct an EFT in which the dark mat- 
ter interacts directly with pairs of electroweak bosons. 
Given our assumption that x is a SM gauge singlet, all 
such interactions are higher dimensional operators. Such 
operators begin at dimension 7, though through elec- 
troweak symmetry breaking they also imply effectively 
dimension 5 descendant operators as well. 

The dimension 5 terms originate from. 
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XX 
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where D^H is the ordinary covariant derivative acting 
on the SM Higgs doublet. Expanding out the covariant 
derivative and replacing H by its vacuum expectation 
value, we arrive at 



(3) 



ft is worth noting that while the overall size of both cou- 
plings may be varied by shifting u^/Ag, the ratio of the 
couplings to pairs of W and Z bosons are fixed with re- 
spect to one other. At higher order, this operator also 
results in couplings to pairs of photons and to Zj through 
loops of W bosons. 

At dimension 7, there are also couplings to the kinetic 
terms of the electroweak bosons. 
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FIG. 1: Representative diagrams for production of dark mat- 
ter pairs (xx) associated with a Z boson in theories where 
dark matter interacts with quarks (top) or directly with Z 
boson pairs (bottom). 



TABLE I: Production cross sections (in fb) for pair produc- 
tion of WIMPs in association with a Z boson, pp — > Zxx ~^ 
£"'"£~XXi in theories where the dark matter interacts primarily 
with quarks, for A, = 1 TeV and ^ = 7 TeV. 
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where Sw and c^, are the sine and cosine of the weak mix- 
ing angle, respectively. For these kinetic operators, the 
over- all size can be though of as controlled by k2/Aj, but 
there is still freedom to adjust the relative importance of 
various pairs by adjusting ki/k2. 



where Fi, i = 1,2,3 are the field strengths for the SM 
f/(l), SU{2), and SU{3) gauge groups. The couphngs of 
dark matter to pairs of SM gauge bosons are given by: 
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DARK MATTER PRODUCTION IN 
ASSOCIATION WITH A Z BOSON 



The process of interest is pair-production of dark mat- 
ter particles in conjunction with one Z boson. Represen- 
tative Feynman diagrams are shown in Figure [T] for both 
the case of interactions with quarks as well as dark mat- 
ter which interacts directly with weak bosons. In order 
to match on to the existing ATLAS ZZ measurement, 
we consider collisions at s — 1 TeV. Cross sections 
for each of the quark operators of Eq. ([!]) are presented 
for various dark matter masses and the corresponding A 
set equal to 1 TeV in Table |Ij These cross sections scale 
as cx 1/Aqj^ for operator Dl, and as oc 1/A'* for D5, D8, 
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and D9. The rates for Dl are considerably smaller, due 
to suppression of the contribution of valence quark con- 
tributions by the small up and down quark masses. For 
this operator, loops involving the top quark are expected 
to be sizable and will result in an increase of the rate 

m- 

In part because of the rather stiff branching ratio 
penalty in asking for Z — >■ e+e" or /x^/i", the rate for 
the operators describing interactions with quarks is some- 
what smaller than the corresponding rates for mono-jets, 
mono-photons, or mono-VFs. However, the mono-Z sig- 
nature is nonetheless worth exploring, in part because 
it samples a different weighting of couplings to up-type 
versus down-type quarks, but also because the system- 
atic uncertainties on the backgrounds should scale more 
favorably for mono-Z than for mono-jets or even mono- 
photons, given that fake "QCD" backgrounds should be 
much smaller for mono-Zs. 

For the case of direct interactions with weak bosons, 
the dimension 5 operator is mediated only by Z exchange, 
whereas the dimension 7 operator contains a mixture of 
Z and photon exchange, with the relative importance of 
the two controlled by fci/fc2- We consider two example 
admixtures of the dimension 7 operators: 

• ki — k2 leading to contributions from both Z and 
7 exchange. 

• ki — c^/s^fc2, for which the 7 exchange graph is 
negligible. 



ATLAS RESULTS 

The ATLAS collaboration performed a measurement of 
the ZZ — )■ Uvv cross-section [Sj in pp data with y/s = 7 
TeV and integrated luminosity of 4.6 fb^^. The fiducial 
region is defined as: 

• two same-flavor opposite-sign electrons or muons, 
each with > 20 GeV, |7?^| < 2.5; 

• dilepton invariant mass close to the Z boson mass: 
mu e [76, 106] GeV; 

• no particle-level jet with plj, > 25 GeV and 
W\ <4.5; 

• (br-pfl)/pf <0.4; 

• -p^^ X cos(A0(p^^pf )) > 75 GeV. 

The results are consistent with Standard Model expec- 
tations (within Icr for the n'^ fJ-~ channel and within 2a 
over-all), as shown in Table [TT} 

We use the expected background yield with uncertain- 
ties to calculate an upper limit on the number of events 
due to a new source which could be present in the col- 
lected data. Using the CLs method [551 US]; we find 



TABLE IL Expected backgrounds and observed data in the 
ATLAS ZZ Uvu analysis ^8, in pp collisions at a/s = 7 TeV 
with integrated luminosity of 4.6 fb~^. The first uncertainty 
is statistical and systematic and the second uncertainty is 
luminosity. 
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N < 18.0 at 90% confidence limit (CL). We convert this 
to a limit on the cross section (cr) for new physics based 
on N = a X e X C, where e is the fraction of new physics 
events which satisfy the selection requirements. 

To aid the reinterpretation of the result, ATLAS has 
divided their calculation of e = Azz x Czz into two 
pieces: the fiducial acceptance (Azz), the fraction of 
events which fall into a specified parton-level fiducial re- 
gion, and the reconstruction efficiency {Czz), the frac- 
tion of events in the fiducial region which satisfy the final 
selection. We expect the reconstruction efficiency Czz 
to be largely model independent, as the fiducial region is 
chosen such that Czz is determined by the detector per- 
formance for specific final state objects rather than the 
production mechanism. Therefore, a calculation of the 
fiducial acceptance for a new model is all that is needed 
for an estimate of the total efficiency e. In terms of the 
fiducial acceptance Azz, the limits are 



a < 



N 



Azz X Czz X C 
which results at the 90% CL in, 

ICQ 

(7(90% CL) < 



Azz X 0.679 X 4.6 fb~^ 



(10) 



(11) 



INTERPRETATION 

We simulate the predicted Zxx events for each of the 
EFTs described above, using MADGRAPH [27], with show- 
ering and hadronization provided by pythia [5H] and 
particle-level jet clustering with fast jet [^Hj- We work 
at tree level, though it should be noted that very recently 
the next to leading order rates for several operators have 
been computed [30], and result in a modest increase in 
the expected rates. For each of the EFTs and a variety of 
dark matter masses, we compute the fiducial acceptance 
as defined above. The resulting acceptances for each EFT 
as a function of the dark matter mass are shown in Fig- 
ure [2| The acceptance varies between about « 20 — 35%, 
depending on the operator, and is roughly constant up 
to dark matter masses of about 1 TeV. 

Combined with the ATLAS measurement, the accep- 
tances translate into bounds on the production cross sec- 
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FIG. 2: Acceptance of the dark matter production process, 
pp Zxx, with Z e+r, for the ATLAS ZZ Uw 
fiducial region (see text) as a function of the dark matter 
mass. 



FIG. 3: Distribution of axial missing transverse momentum 
in simulated Standard Model ZZ — >■ Uvu and pp — >■ Zxx^ in 
the EFTs described in the text, for A = 1 TeV and = 100 
GeV at = 7 TeV. 



tion of a{Zxx) ^ 10 - 100 fb at 90% CL. It is worth 
noting that the fiducial acceptance is significantly higher 
for dark matter production than for the SM ZZ produc- 
tion due to the larger missing transverse momentum in 
the ZxXi ^ shown in Figure [3] for a sample parameter 
point with = 100 GeV and A = 1 TeV. For each 
interaction type, we determine the lower bound on the 
scale A which characterizes its strength. The results are 
presented in Figure |4j which indicate that for some types 
of interactions, scales on the order of 100 GeV to TeV 
can be probed. 

For the case of interactions with quarks, we further 
translate these bounds into the plane of spin- independent 
(SI) and spin-dependent (SD) scattering with nucleons, 
as related in [TH| . These bounds are shown in Figure [s] 
along with bounds from CoGeNT [3T] and Xenon 100 
[32j . and competing collider results from mono-jet and 
mono-photon searches. As is typical, bounds from col- 
liders provide a unique probe of very light dark matter 
particles, and dominate as probes of spin-dependent in- 
teractions. Of course, the collider bounds are subject to 
the assumption that the EFT containing a contact inter- 



action is a good description of the physics. In cases with 
light mediating particles, these bounds can sometimes be 
weaker. Over-all, the two searches exhibit a high degree 
of complementarity. 

CONCLUSIONS 

In this article, we have looked at the collider produc- 
tion of dark matter pairs in association with a Z boson 
which decays into charged leptons, a signature we refer 
to as a mono-Z. We work in the context of effective field 
theories in which the dark matter either interacts directly 
with quarks, or with a pair of electroweak bosons. We 
derive limits on the strength of such interactions based 
on the recent ATLAS measurement of ZZ production 
(where one Z decays into charged leptons and the other 
into neutrinos) ands find that the current limits already 
probe the TeV scale for some types of interactions. 

For the case of interactions directly with quarks, the 
mono-Z signature provides limits which are somewhat 
weaker than those from mono-jets or mono-photons. 
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FIG. 4: Excluded region (blue) at 90% CL for the indicated 
interactions as a function of dark matter mass. The two re- 
gions for the dimension 7 ZZxx correspond to the choices of 
k\/k2 discussed in the text, with either maximal contribution 
from photon graphs (upper curve) or no contribution (lower 
curve) . 



FIG. 5: Exclusion regions at 90% CL in the plane of the 
dark matter mass (m^) versus dark matter - nucleon spin- 
independent (top) or spin-dependent (bottom) cross-section 
plane. Results from this analysis are compared to existing 
collider and direct detection limits. 



Nonetheless, mono-Z searches are expected to be less 
subject to systematic uncertainties from jet energy scales 
and photon identification, and thus may scale better at 
large luminosities. If a discovery is made, the mono-Z 
signature offers a different way to dissect the couplings 
of up-type versus down-type quarks. If the dominant in- 
teraction is instead to pairs of weak bosons, colliders offer 
a unique opportunity for discovery. 

Our results illustrate the complementarity between col- 
lider and direct searches of dark matter, and show how 
together they result in a more complete picture of dark 
matter interactions with the SM fields. 
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